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g ■ Abstract 

■ We give a complete geometrical description of the effective Hamiltonians common in nuclear shell 

model calculations. By recasting the theory in a mamfestly geometric form, we remterpret and 
\ clarify several points. Some of these results are hitherto unknown or unpublished. In particular, 

commuting observables and symmetries are discussed in detail. Simple and explicit proofs are 
given, and numerical algorithms are proposed, that improve and stabilize common methods used 
today. 
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I. INTRODUCTION 



Effective Hamiltonians and interactions are routinely used in shell-model calculations of 
nuclear spectra . The published mathematical theory of the effective Hamiltonian 



IS com 
etc 



implicated and usually focuses on perturbation the 
• fl 3 y, 3, y, y, [W, in, 12, lal. in t^s artier, 



icated and usually focuses on perturbation theoretical aspects, diagram expansions, 

we recast and reinterpret the basic 
elements of the theory geometrically. We focus on the geometric relationship between the 
exact eigenvectors \il)k) and the effective eigenvectors both for the usual non-Hermitian 

and for the Hermitian effective Hamilto- 



Bloch-Brandow effective Hamiltonian 
nian 
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15j, which we dub the canonical effective Hamiltonian due to its geometric 



significance. This results in a clear geometric understanding of the de-coupling operator 
uj = QujP, and a simple proof and characterization of the Hermitian effective Hamilto- 
nian in terms of subspace rotations, in the same way as the non-Hermitian Hamiltonian is 
characterized by subspace projections. 

As a by-product, we obtain a simple and stable numerical algorithm to compute the exact 
effective Hamiltonian. 

The goal of effective interaction theory is to devise a Hamiltonian H^s in a model space 
V of (much) smaller dimension m than the dimension n of Hilbert space Ti, with m exact 
eigenvalues of the original Hamiltonian H = Hq + Hi, where Hi is usually considered as a 
perturbation. The model space V is usually taken as the span of a few eigenvectors {|efc)}^;^ 
of Hq, i.e., the unperturbed Hamiltonian in a perturbational view. 

Effective Hamiltonians in A-hodj systems must invariably be approximated (otherwise 
there would be no need for Hes), usually by perturbation theory, but a sub-cluster approxi- 
mation is also possible 3|, |6| . In that case, the exact a-body canonical effective Hamiltonian 
is computed, where a < A. From this, one extracts an effective a-body interaction and apply 
it to the y4-body system. In this case, we present a new algorithm for computing the exact 
effective interaction that is conceptually and computationally simpler than the usual one 
which relies on both matrix inversion and square root 3|, lUl] , as the only non-trivial matrix 
operation is the singular value decomposition (SVD). 

The article is organized as follows. In Sec. [Ill we introduce some notation and define the 
singular value decomposition of linear operators and the principal angles and vectors be- 
tween two linear spaces. In Sec. Illll we define and analyze the Bloch-Brandow and canonical 
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effective Hamiltonians. The main part consists of a geometric analysis of the exact eigen- 
vectors, and forms the basis for the analysis of the effective Hamiltonians. We also discuss 
the impact of symmetries of the Hamiltonian, i.e., conservation laws. In Sec. IIVI we give 
concrete matrix expressions and algorithms for computing the effective Hamiltonians, and 
in the canonical case it is, to the author's knowledge, previously unknown. In Sec. |V] we 
sum up and briefly discuss the results and possible future projects. 

II. TOOLS AND NOTATION 
A. Linear spaces and operators 

We shall use the Dirac notation for vectors, inner products and operators, in order to 
make a clear, basis-independent formulation. By JF, etc., we denote (finite dimensional) 
Hilbert spaces, and vectors are denoted by kets, e.g., {ip), as usual. Our underlying Hilbert 
space is denoted by Ti, with n = dim(7i). In general, n is infinite. We shall, however, assume 
it to be finite. Our results are still valid in the infinite dimensional case if H is assumed to 
have a discrete spectrum and at least m linearly independent eigenvectors. 

We are also given a Hamiltonian H, a linear, Hermitian operator (i.e., H = if^) on Ti. 
Its spectral decomposition is defined to be 

n 

H = J2Ek\i^k){i^k\. 

k=l 

Thus, Ek and \ipk) are the (real) eigenvalues and (orthonormal) eigenvectors, respectively. 

We are also given a subspace V C H, called the model space, which in principle is 
arbitrary. Let {|efc)}'^]^ be an orthonormal basis, for definiteness, viz, 

P := span{|eA;) : k = l,---,m}. 

Let P be its orthogonal projector, i.e., 

m 

P-.n^V, P = ^|ej)(ej|, m = dim(P)<n, 
i=i 

The basis {\ej)}'jLi is commonly taken to be eigenvectors for Ho. 

The orthogonal complement of the model space, Q = V^, has the orthogonal projector 
Q = I — P, and is called the excluded space. 
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This division of 7i into V and Q transfers to operators in Ti. These are in a natural way 
spht into four parts, viz, for an arbitrary operator A, 

A = {P + Q)A{P + Q)= PAP + PAQ + QAP + QAQ, (1) 

where PAP maps the model space into itself, QAP maps V into Q, and so forth. It is 
convenient to picture this in a block-form of A, viz, 

^ PAP PAQ 
QAP QAQ 

B. The singular value decomposition 

A recurrent tool in this work is the singular value decomposition (SVD) of an operator 
A : X —>■ y. Here, p = dim{X) and q = dim(3^) are arbitrary. Then there exists orthonormal 
bases {\xk)}^^i and {\yk)}l=i of X and 3^, respectively, and r = min(p, g) non-negative real 
numbers with ak > (Jk+i for all k, such that 

r 

A = '^crk\yk){xk\. 

k=l 

This is the singular value decomposition (SVD) of A, and it always exists. It may happen 
that some of the basis vectors do not participate in the sum; either if p ^ q, or if cx/c = for 
some k. 

The vectors \xk) are called right singular vectors, while \yk) are called left singular vectors. 
The values are called singular values, and A is one-to-one and onto (i.e., nonsingular) if 
and only if cxfc > for all k, and p = q. The inverse is then 




as easily verified. 

A recursive variational characterization of the singular values and vectors is the following 
Q: 



ak = max max Re(f|yl|M) 

\u)ex, {u\u)=i \v)ey, {v\v)=i 

{u\uj)=0,j<k {v\vj)=0,j<k 

= ■■ {vk\A\uk). (2) 
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The latter equality implicitly states that the maximum is actually real. The SVD is very 
powerful, as it gives an interpretation and representation of any linear operator A as a 
simple scaling with respect to one orthonormal basis, and then transformation to another. 
The singular vectors are not unique, but the singular values are. 



C. Principal angles and vectors 

Important tools for comparing linear subspaces T and Q oi Ti are the principal angles 



and principal vectors 
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18j. The principal angles generalize the notion of angles between 



vectors to subspaces in a natural way. They are also called canonical angles. Assume that 

p = dim(jF) > q = dim(^) > 1. 

(If p < q, we simply exchange JF and Q.) Then, q principal angles 9k G [0,7r/2], with 
9k < 9k+i for all k, and the left and right principal vectors \^k) ^ ^ and \rik) G Q are defined 
recursively through 

cos6k = max max Re(f|?7) 

(CI?j>=0,J<fc {v\Vj)=0,j<k 

=■■ i^klVk). (3) 

Again, the last equality implicitly states that the maximum actually is real. One sees that 
9k is the angle between \^k) G and \rik) G Q. 

It is evident from Eqns. ([2]) and ([3]) that the principal angles and vectors are closely 
related to the SVD. Indeed, if we consider the product of the orthogonal projectors Pjr and 
Pg and compute the SVD, we obtain 

p q q 

PrPg = \^k){^k\ \r]k){Vk\ = cos^fc|^fc)(%|, 

k=l j=l k=l 

where we extended the orthonormal vectors {\^k)}'k=i with p — q vectors into a basis for JF, 
which is always possible. This equation in particular implies the additional orthogonality 
relation {^j\rik) = 6j^k cos 9k on the principal vectors. 

The principal vectors constitute orthonormal bases that should be rotated into each other 
if the spaces were to be aligned. Moreover, the rotations are by the smallest angles possible. 
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III. EFFECTIVE HAMILTONIANS 



A. Similarity transforms 

The goal of the effective Hamiltonian is to reproduce exactly m of the eigenvalues, and 
(necessarily) approximately m of the eigenvectors. We shall assume that the first m eigen- 
pairs {Ek, I'lpk)), k = 1, . . . ,m, defines these. We define the space £ as 

£ := span{|'?/'fc) : k = 1, - ■ ■ , m}. 

The orthogonal projector P' onto £ is 

m 

p' = J2m{^pk\. (4) 

k=l 

We denote by {Ek, li'^)), k = 1, . . . ,m, the effective Hamiltonian eigenvalues and eigen- 
vectors. Of course, the G V must constitute a basis for V, but not necessary an 
orthonormal basis. Geometrically, we want \ipl^) to be as close as possible to \ipk), i-e., we 
want £ to be as close to V as possible. 

Let iV'fc^) be the bi-orthogonal basis, i.e., {'ipj^l'ipl^) = 5j,fc, so that 



k=l k=l 

The spectral decomposition of H^s becomes 

m 

H,^ = Y,Ek\rk){t 



k=l 



Since Hes is to have eigenvalues identical to m of those of H, and since Hes operates only 
in V, we may relate H^s to H through a similarity transform, viz, 

H^s = PHP = P{e-^He^)P, (5) 

where exp(S') exp(— S*) = /. Any invertible operator has a logarithm, so Eqn. ([5]) is com- 
pletely general. 

Now, Hcs = PHP is an effective Hamiltonian only if the Bloch equation 

QHP = Qe-^He^P = (6) 
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FIG. 1: Plane spanned by \^k) and \r]k)- Action of projectors P and Q on \rjk) indicated 



is satisfied js], since V is then invariant under the action of H. The eigenvectors of ifes are 
now given by 

l^f ) = e-^l^fc) G P, fc = l,---,m. (7) 

Thus, an effective Hamiltonian can now be defined for every 5* such that Eqn. ([6]) holds. It 
is readily seen that H = W if and only if 5* is skew- Hermit ian, i.e., that 5*^ = —S. There is 
still much freedom in the choice of exponent S. Indeed, given any invertible operator A in 
"P, A~^HciiA is a new effective Hamiltonian with the same effective eigenvalues as ifefr, and 
= A-iexp(-^)lV). 



B. Geometry of the model space 

We will benefit from a detailed discussion of the spaces £ and V before we discuss the 
Bloch-Brandow and the canonical effective Hamiltonians in detail. 

Since dim(P) = dim(£') = m, the closeness of the effective and exact eigenvectors can be 
characterized and measured by the orientation of £ relative to V in 7i, using m canonical 
angles Ok and principal vectors \rjk) G £ and \^k) ^ 'P- Recall, that co?,6k = {C,k\Vk) and that 
the angles 9k G [0,7r/2] were the smallest possible such that the principal vectors are the 
orthonormal bases of V and £ that are closest to each other. 

We now define the unitary operator Z = exp(G') that rotates V into £ according to this 
description, i.e., we should have Z\C,k) = \Vk)- In Fig- [Hthe plane spanned by \rik) and \C,k) 
if 9k > is depicted. Recall, that = cos 9 j6j^k- Note that 1^^) = \rik) if and only if 

9k = 0, and the plane degenerates into a line. If 9^ > 0, the vector \xk) is defined so that it 
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together with \^k) is an orthonormal basis for the plane, viz, 

\Vk) = P\Vk) + Q\Vk) = cos{9k)\^k) + sin(6'fc)|xfc), (8) 

where 

Thus, {Ixfc)} U {l^fc)} is an orthonormal basis for V ®S, whose dimension is 2m — n^, where 
Uz is the number of 6k = 0. The set {\xk) '■ Ok > 0, k = 1, - ■ ■ ,m} is an orthonormal basis 
for QS which contains Q\ipk) for all = 1, ■ ■ ■ , m. 

The operator Z is now defined as a rotation in "P © i.e., by elementary trigonometry, 

Z\^k) ■■= \Vk) 

Z\7]k) := 2cosek\Vk) - \Ck) (10) 
In terms of the orthonormal basis, we obtain a manifest planar rotation for each k, i.e., 

Z\xk) = cos 9 k\xk) - sine k\^k) 

Z\^k) = sin 9k\xk) + cos 9k\^k)- (11) 

On the rest of the Hilbert space, TiQiV ®£)^ Z is the identity. The operator Z implements 
the so-called direct rotation of V into 8. From Eqn. ffTTj) we obtain 



Z := I + Y}^os{9k)-mXk){Xk\ + \^k){ik\) 

k=l 
m 

+ J2M0k){\Xk){^k\-\^k){Xk\), (12) 

k=l 

It is instructive to exhibit the Lie algebra element G G su{n) such that Z = exp(G) G 
SU{n). Since we have Eqn. (|TT|) . is is easy to do this. Indeed, taking the exponential of 

m 

G = J2^k{\Xk){^k\-\^k){Xk\), (13) 
fe=i 

by summing the series for sin(^) and cos(6'), we readily obtain Z = exp{G), the desired 
result. Moreover, observe that the fc'th term in Eqn. f[T^ commutes with the j'th term, so, 
exp(G') is exhibited as a sequence of commuting rotations using the canonical angles 9k- 
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C. The Bloch-Brandow effective Hamiltonian and the decoupHng operator 



For reference, we review some 



which we denote by Hf^ 



_)roperties of the Bloch-Brandow effective Hamiltonian, 
The effective eigenvectors are defined by 



(14) 



Since \Pipk) are the orthogonal projections of onto V, we deduce that the Bloch-Brandow 
effective eigenvectors are the closest possible to the exact model space eigenvectors. In this 
sense, the Bloch-Brandow effective Hamiltonian is the optimal choice. 

It is obvious that is non-Hermitian, as rejecting the excluded space eigenvector 

components renders the effective eigenvectors non-orthonormal, i.e., 

(^f|V^f) = 5,-,-(V^,|g|^,)^5,-,. 



In terms of similarity transforms, we o 
operator or correlation operator [ll] 



3tain HffP by setting S = u, the so-called de-coupling 
It is defined by u; = QuP and the equation 



iuPliJk) ■■=Q\^k). 

Again, for this to be a meaningful definition, {\Pipk)}k=i niust be a basis for V. 
Since tu^ = 0, exp(±u;) = 1 ± u;, and Eqn. (JTj) becomes 

e-^\^k) = {l-oo)\^k) = (1 - QMk) = IPi^k). 



(15) 



For if™ we thus obtain 



H^^"" = Pe-'^He'^P = PH{P + lj) 



(16) 



After this initial review, we now relate u to the geometry of S and V. The SVD of 
uj is readily obtainable by expanding the principal vectors {\r]k)}^=i in the m eigenvectors 
{\ipk)}k=i^ sets which both constitute a basis for S, and inserting in Eqn. (fT5l) . We have 

m 

Q\^k) = ^Q\'^j){^j\rik) 

m 
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that is, 

Lo {co?,9k\ik)) = sin6'fc|Xi)- 

The result is 

m 

= ^tan6'fc|xfc)(6|, (17) 

k=l 

which is the SVD of u. The operator u is thus exhibited as an operator intimately related 
to the principal angles and vectors of V and S: It transforms the principal vectors of V into 
an orthonormal basis for QS, with coefficients determined by the canonical angles 6k- Using 
Eqn. ([8]) we obtain an alternative expression, viz, 

m ^ 

^cos^, 

D. The canonical effective Hamiltonian 



Hermitian effective Hamiltonians have independently been introduced by various au- 

introduced a unitary transformation H = 



thors since 1929, when Van Vleck 
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exp(-5) 



Primas 



H exp{S) to decouple the model space to second order in the interaction. In 1963, 



2ii 



considered an order by order expansion of this H using the Baker-Campbell- 



Hausdorff formula and commutator functions to determine S, a technique also used in many 
other settings in which a transformation is in a Lie group, see, e.g.^Ref. {22!] and refer- 
ences therein. This approach was elaborated by Shavitt and Redmon [7i], who were the first 
to mathematically connect this Hermitian effective Hamiltonian to H^^, as in Eqn. ( l27|l 
below. In the nuclear physics community, Suzuki [23] has been a strong advocate of Hermi- 
tian effective interactions and the a-body sub-cluster approximation to the A-hodj effective 



interaction |3|, [111, |23| . Hermiticity in this case is essential. 

Even though a Hermitian effective Hamiltonian is not unique due to the non-uniqueness 
of S* = —S'^, the various Hermitian effective Hamiltonians put forward in the literature all 
turn out to be equivalent [6|]. In the spirit of Klein and Shavitt we employ the term 

"canonical effective Hamiltonian" since this emphasizes the "natural" and geometric nature 
of the Hermitian effective Hamiltonian, which we denote by H^f^. 
Recall the spectral decomposition 



eflfl 
k h 



k=l 
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where the (orthonormal) effective eigenvectors are now defined by the following optimization 
property: The effective eigenvectors {ipl^) are the closest possible to the exact eigenvectors 
\ipk) while still being orthonormal. Thus, where the Bloch-Brandow approach globally mini- 
mizes the distance between the eigenvectors, at the cost of non-orthonormality, the canonical 
approach has the unitarity constraint on the similarity transformation, rendering H^^ Her- 
mitian. 

Given a collection $ = {|</)i), C V oi m vectors, which are candidates for 
effective eigenvectors, define the functional S'[$] by 

m 

sm := Elll'^'^)-!^'^)!!' 

k=l 

m 

= K+j2\mk)-\M\' 

k=l 

m m 

= K + J2\\\<Pk)f - 2ReJ2{^k\P\<Pk), (19) 

k=l k=l 

where K = YlT=i \\Q\i'k)\\'^ is a constant in this context. The effective eigenvectors are now 
minimizers of S[^]. 

The global minimum, when $ C P is allowed to vary freely, is attained for \(j)k) = {P'ipk), 
the Bloch-Brandow effective eigenvectors. However, the canonical effective eigenvectors are 
determined by minimizing 5'[$] over all orthonormal sets $, which then becomes equivalent 
to maximizing the last term in Eqn. (1191) . i.e., the overlaps Re('?/'fc|P|0fc) under the 
orthonormality constraint. 

We will now prove the striking fact that the solution is given by 

lV'f) = l0.) = e-^|^/^,), (20) 

where the unitary operator Z := exp(G') e SU{n) is the rotation (fT2|) . Equation (120|) should 
be compared with Eqn. ([7]). Thus, the exact eigenvectors are simply the direct rotations of 
the effective eigenvectors from the model space into S. 

Let us expand lipk) G ^ and G "P in the principal vector bases, viz, 

m 

l^fc) = ^\vj){vj\i^k), 

i=i 

m 

11 



Using {rjjl^k) = Sj^kCOsOj, we compute the sum A := ^'k^=i{i^k\P\(pk) as 

m 

^ = J2 {MVj){Vj\Ci){Ce\<Pk) 

k,j,£=l 
m 

= cosej{^j\(t)k){ipk\vj) 
j,k=i 

Now, 

m 

where Uj^k is a unitary matrix, which imphes \ujj\ < 1. Moreover, Ujj = 1 for all j if and 
only if Uj^k = ^j,k, which then maximizes A, and also lie A. Thus, 

m 

Y{^j\<Pk){ipk\ve) = Sj,e, 

k=l 

i.e., 

{^j\<Pk) = (Vjl-ipk) = {^j\Z^ipk), 

from which Eqn. (!20|) follows since (1^^)}^]^ is a basis for V, and the proof is complete. 
The similarity transform in Eqn. is thus manifest, with S = G, viz, 

H^s = PZ^HZP = Pe-^He^P. (21) 

Moreover, QHP = PHQ = 0, verifying that the direct rotation in fact block diagonalizes 
H. 

E. Computing \ipf.^) 

Assume that \Pipk) '■= Pli^k), k = 1, ■ ■ ■ ,m are available. The effective eigenvectors 
are then given by a basis change F, i.e., the operator F : V V defined by 

F|P^,) := l^f ). 
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Using the principal vector basis we obtain 

m 

rn 

= COS ej\Q{r]j\iJk) 

m 
m 

from which we get the SVD 




= {uj^uj + Py/', (22) 

where we have used Eqn. (ITSl) . From Eqn. (!22l) we see that F is symmetric and positive 
definite. Moreover, smaller angles 6k means F is closer to the identity, consistent with S is 
closer to V. 

Let {Pipk) now be given in the orthonormal "zero order" basis {|efc)}^^ for V, i.e., we 
have the basis change operator U given by 

m 

U:=Y,\P4^k){ek\, (23) 
fc=i 

which transforms from the given basis to the Bloch-Brandow effective eigenvectors. In terms 
of the principal vector basis, 

m 

U = ^cos9j\Q{r]j\^\^k){ek\ 

j=l k 
m 

= : J2^ose,\Q{y,l (24) 
i=i 

which is, in fact, the SVD since the last sum over is a unitary map from V to S. In the 
operator UW this basis-dependent factor cancels, viz, 

m 

k=l 
m 

= ^cos^6'fc|^fc)(^fc|, 

k=l 
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that is, 

F = {UU^)-^'\ 

If we seek \'ip'^) in the basis {|efc)}^i as well, we let V be the corresponding basis change 
operator, i.e., 

V ■=FU = [UU^Y^I'^tj. (25) 

Equation (125|) shows that \ip'^) is obtained by "straightening out" l-PV'fc); and that this 
depends only on the latter vectors. This is, in fact, an alternative to the common Gram- 
Schmidt orthogonalization used in mathematical constructions and proofs. It was first intro- 



duced by Lowdin 



24| under the name "symmetric orthogonalization" , and so-called "Lowdin 



bases" are widely-used in quantum chemistry, where non-orthogonal basis functions are or- 
thogonalized according to Eqn. ( l25i) . It seemingly requires both inversion and matrix square 
root, but is easily computed using the SVD. Combining Eqns. fl22|) and fl24|) gives 

m 

V = Y.\^^)^y^\^ (26) 
fc=i 

so that if the SVD ( 12^ is available, V is readily computed. Eqn. (l26l) is easily expressed in 
terms of matrices, but we defer the discussion to Sec. HVl 



F. Shavitt's expression for exp(G) 

Shavitt and Redmon [?! proved that 

G = tanh-i(cj-cut) (27) 

gives the Lie algebra element for the unitary operator Z = exp(G). The quite complicated 
proof was done using an expansion of the similarity transform using the Baker-Campbell- 
Hausdorff formula. 

It may be clear now, that in the present context we obtain the result simply as a by- 
product of the treatment in Section IIII Bl and the SVD 0171) of u, given in terms of the 
principal vectors and angles. We prove this here. 

The function tanh~^(z) is defined by its (complex) Taylor expansion about the origin, 
i.e., 

tanh-H^) = 5^ ^--y. (28) 



2n + l 

n=0 
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The series converges for \z\ < 1. Moreover, 



2 \1- z 

also valid for \z\ < 1. For z := uj — uj'' we compute 



tanh-^(z) = ilnfi±^K (29) 



k=l 

Using orthogonality relations between \C,k) and \xk) we obtain 

m 

^2n+l ^ i-irY,f^l''^\\Xk){^k\-\^k){Xk\) 
k=l 

m 

= ^Y.H^ikr^\\xk){^k\-\^k){xk\). 

k=l 

Using i tanh~"^(— Z2;) = tan~^(z), we sum the series (l28l) to 

m 

G = tanh-i(z) = Y^9k{\xk){^k\ - \^k){Xk\), 

k=l 

which is identical to Eqn. (fT3!) . The series does not converge for 6k > vr/4, but the result is 
trivially analytically continued to arbitrary < 6k < 7^/2. 

Ill ] to compute ifgg in 



We now turn to the effective Hamiltonian. It is common 
terms of u directly, using the definition ( l29i) of tanh~^(z), which implies 



g±tanh-i(2) II ±Z 1±Z 



Upon insertion into H = exp{—G)Hexp{G), we obtain 

~ 1 — UJ + Uj'^ 1 + u — Uj"^ 

H = — — H - 



Projecting onto V, the effective Hamiltonian becomes 

H,s ={P + uo^uo)-^'^{P + uo^)H{P + cu)(P + uo^uoY^/^. (30) 

By using the Bloch-equation ([6]) for the Bloch-Brandow effective Hamiltonian, we may elim- 
inate QHQ from the above expression for iJcfr? yielding 

His = (P + u^ujf'^H{P + uj){P + uj^uj)-"\ (31) 
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This expression is commonly implemented in numerical applications [3|, |25| . By comparing 
with Eqns. f l22l) and f|T6|) we immediately see that 

H:^ = FH^iF-\ (32) 

which gives H^^ as a similarity transform of . In themselves, Eqns. fl3Tl) and fl32l) are not 
manifestly Hermitian, stemming from the elimination of QHQ. An implementation would 
require complicated matrix manipulations, including a matrix square root. It is therefore 
better to compute using 

k 

together with Eqn. (126!) . where the most complicated operation is the SVD of the operator 
U given by Eqn. fl23|l . In Sec. [IV]we give a concrete matrix expression for Hl^. 



G. Commuting observables 

Great simplifications arise in the general quantum problem if continuous symmetries of 
the Hamiltonian can be identified, i.e., if one can find one or more observables S such that 
[if, S] = 0. Here, we discuss the impact of such symmetries of H on the effective Hamiltonian 
Hes] both in the Bloch-Brandow and the canonical case. We point out the importance of 
choosing a model space that is an invariant of S as well, i.e., [S,P] = 0. In fact, we prove 
that this is the case if and only if [ifeff, "S*] = 0, i.e., ifes has the same continuous symmetry. 

Let = S*^ be an observable such that [H, S] = 0, i.e., H and S have a common basis of 
eigenvectors. We shall assume that {\i/jk)}k=i is such a basis, viz, 

S\^k) = Skl^k)- (33) 

In general, there will be degeneracies in both and Sk- 
We now make the important assumption that 

[S, P] = 0, (34) 

which is equivalent to 

5* = PSP + QSQ. 
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Under the assumption we have 

S\Pi;k)=PS\iJu)=Sk\Pi^u), 

so that the Bloch-Brandow effective eigenvectors are still eigenvectors of S with the same 
eigenvalue s^- Moreover, as we assume that {\Pil)k)}^=i is a (non-orthonormal) basis for P, 
this not possible if [S, P] ^ 0. Thus, [H^i , ^] = if and only if [S, P] = (in addition to 
the assumption (!33|) .) 

The assumption flMl) also implies that [5, cu] = 0, where oj = QujP is the de-coupling 
operator. We prove this by checking that it holds for all \il)k)- For k < m, 

ujS\ipk) = SkQlipk) = SQlipk) = Suj\ipk), (35) 
while foT k > m we need to expand Plipk) in \P'^j), j ^ ""^j viz, 

m 

P\iJk) = \P^j){Pi^j\P\^k), k>m 

and use Eqn. (!35|) . Furthermore, [5', cj''']"!' = [c<j''",S'] = 0. It follows that 

[S,{uj -uj^Y] = Q, n = 0,l,..., 

and, by Eqn. ([2S]), that 

[S, e«] = [5, e-^] = 0. 

This gives 

Again, since {|'i/'fc'^)}fcLi is a basis for V, this holds if and only if [5, P] = 0. Accordingly, 
[HIq, S*] = if and only if [S, P] = (and the assumption (15^ .) 

The importance of this fact is obvious. If one starts with a Hamiltonian that conserves, 
say, angular momentum, and computes the effective interaction using a model space that is 
not an invariant for the angular momentum operator, i.e., not rotationally symmetric, then 
the final Hamiltonian will not have angular momentum as a good quantum number. 

One possible remedy if [P, S] 7^ is to define the effective observable S^e ■= 
P exp{—G)S exp{G)P (which in the commuting case is equal to PSP) which obviously 
commutes with Hcs and satisfies 

S,s\^f) = s,\^pf). 
17 



This amounts to modifying the concept of rotational symmetry in the above example. 

The assumptions fl33l) and fl34|) have consequences also for the structure of the principal 
vectors ^ V and \rik) G S. Indeed, write 

S = 

s 
s 

where the sum runs over all distinct eigenvalues Sk, k = 1, ■■■ ,m of S, and where Ss (Vs) is 
the corresponding eigenspace, i.e., 

£s := span{\ipk) ■ S\ipk) = slipk)} 
Vs := spanjl^f) : S\ijf) = s\i^f)}. 

The eigenspaces are all mutually orthogonal, viz, Sg-^^s', Vs-i-Vs', and Ss-^Vs', for s ^ s' . 
The definition ([3]) of the principal vectors and angles can then we written 

cos(^fc) = max max max Re(^|r7) 

^ lOePs, (C|?)=l \ri)(i£s,{v\r})=l 
felO=0,i<fc (r,j\ri)=0,j<k 

Thus, for each k, there is an eigenvalue s of such that 

S\Vk) = s\r]k), 

showing that the principal vectors are eigenvectors of S if and only if [S, P] = 0, [S, H] = 0, 
and the assumption fl33|) . 

The present symmetry considerations imply that model spaces obeying as many symme- 
tries as possible should be favored over less symmetric model spaces, since these other model 
spaces become less "natural" or "less effective" in the sense that their geometry is less similar 
to the original Hilbert space. This is most easily seen from the fact that principal vectors 
are eigenvectors for the conserved observable 5*. This may well have great consequences 
for the widely-used sub-cluster approximation to the effective Hamiltonian in no-core shell 
model calculations , where one constructs the effective Hamiltonian for a system of 

a particles in order to obtain an approximation to the A > a-body effective Hamiltonian. 

18 



The model space in this case is constructed in different ways in different implementations. 
Some of these model spaces may therefore be better than others due to different symmetry 
properties. 

IV. MATRIX FORMULATIONS 
A. Preliminaries 

Since computer calculations are invariably done using matrices for operators, we here 
present matrix expressions for H^f^ and compare them to those usually programmed in the 
literature, as well as expressions for Hf^ and u. 

Recall the standard basis {|efc)}^^^ of H, where the {|efc)}^^ constitute a basis for V. 
These are usually eigenvectors of the unperturbed "zero order" Hamiltonian Hq, but we will 
not use this assumption. As previously we also assume without loss that the eigenpairs we 
wish to approximate in H^s are {\ilJk)}T=i- 

An operator A -.Ti ^ Ti has a matrix A G C"^" associated with it. The matrix elements 
are given by IKjk = {ej\A\ek) such that 

n n 

^ = XI |ei)(ej|A|efc)(efc| = ^ \ej)^jk{ek\■ 

j,k=l j,k=l 

Similarly, any vector |0) G has a column vector G C" associated with it, with (pj = 
{ej\(f)). We will also view dual vectors, e.g., as row vectors. 

The model space V and the excluded space Q are conveniently identified with C'" and 
C"^™, respectively. Also note that PAP, PAQ, etc., are identified with the upper left 
mxm, upper right m x [n~m), etc., blocks of A as in Eqn. ([1]). We use a notation inspired 
by FORTRAN and matlab and write 

PAP = A(l:m, l:m), PAQ = A(l:m,m + l:n), 

and so forth. 

We introduce the unitary operator U as 

n 

u = Y,m{ek\, 

k=l 
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i.e., a basis change from the chosen standard basis to the eigenvector basis. The columns of 
U are the eigenvectors' components in the standard basis, i.e., 

and are typically the eigenvectors returned from a computer implementation of the spectral 
decomposition, viz, 

H = UEU^ E = diag(Ei,...,E„). (36) 

The SVD is similarly transformed to matrix form. The SVD defined in Sec. Ill Bl is then 
formulated as follows: For any matrix A G C^^^ there exist matrices X G C^^^ {p = min(g, r)) 
and Y G C^'^^, such that X^X = Y^Y = Ip (the identity matrix C^^f), and a non-negative 
diagonal matrix E G M^^^ such that 

A = XEY^ 

Here, T. = diag(cri, ■ ■ ■ , cTp), cjfc being the singular values. 

The columns of X are the left singular vectors' components, i.e., Xj^k = {^j\xk), and simi- 
larly for Y and the right singular vectors. The difference between the two SVD formulations 
is then purely geometric, as the matrix formulation favorizes the standard bases in X and 

The present version of the matrix SVD is often referred to as the "economic" SVD, since 
the matrices X and Y may be extended to unitary matrices over and C, respectively, by 
adding singular values ak = 0, k > m. The matrix Z is then a g x r matrix with "diagonal" 
given by 0"^. This is the "full" SVD, equivalent to our basis- free definition. 

B. Algorithms 

Let the m eigenvectors {ipk) be calculated and arranged in a matrix U, i.e., i/^k = U(l : n, k) 
(where the subscript does not pick a single component). Consider the operator U defined 
in Eqn. fl23l) . whose matrix' columns are the Bloch-Brandow effective eigenvectors \Pipk) in 
the standard basis, viz, 

U = U(l:m, l:m). 

The columns of the matrix oi V = {tJU'^)^^l'^lJ are the canonical effective eigenvectors 'if'^. 
The SVD ([24D can be written 

= xzY^ 
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which gives 

Since Tkk = cos 6k > 0, we obtain 

(0Gt)-i/2 ^xz-^xt, 

which gives, when apphed to U 

V = (uut)-^/20 = XY^ 

Thus, we obtain the canonical effective eigenvectors by talcing the matrix SVD of U = U(l: 
m,l:m) and multiplying together the matrices of singular vectors. As efficient and robust 
SVD implementations are almost universally available, e.g., in the lapack library, this 
makes the canonical effective interaction much easier to compute compared to Eqn. ( l3Ti) . 
viz, 

H;;^ = VE(l:m,l:m)V^ 

This version requires one SVD computation and three matrix multiplications, all with mxm 
matrices, one of which is diagonal. Equation (13T|) requires, on the other hand, several more 
matrix multiplications, inversions and the square root computation. The Bloch-Brandow 
effective Hamiltonian is simply calculated by 

Hfff = OE(l:m,l:m)U-\ 

For the record, the matrix of u is given by 

uj = U(m + l:n, l:m)U^^, 

although we have no use for it when using the SVD based algorithm. It may be useful, 
though, to be able to compute the principal vectors for V and S. For this, one may compute 
the SVD of u or of PP' = UU(1 : n, 1 : mY, the latter which gives cosOk, \C,k) and jr/fc) directly 
in the standard basis as singular values and vectors, respectively. 

V. DISCUSSION AND AND OUTLOOK 

We have characterized the effective Hamiltonians commonly used in nuclear shell-model 
calculations in terms of geometric properties of the spaces V and S. The SVD and the 
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principal angles and vectors were central in the investigation. While the Bloch-Brandow 
effective Hamiltonian is obtained by orthogonally projecting £ onto V, thereby globally 
minimizing the norm-error of the effective eigenvectors, the canonical effective Hamiltonian 
is obtained by rotating £ into V using exp(— G), which minimizes the norm-error while 
retaining orthonormality of the effective eigenvectors. Moreover, we obtained a complete 
description of the de-coupling operator uj in terms of the principal angles and vectors defining 
exp(G). 

An important question is whether the present treatment generalizes to infinite dimen- 
sional Hilbert spaces. Our analysis fits into the general assumptions in the literature, being 
that n = dim(7Y) is large but finite, or at least that the spectrum of H purely discrete. A 
minimal requirement is that H has m eigenvalues, so that £ can be constructed. In partic- 
ular, the SVD generalizes to finite rank operators in the infinite dimensional case, and are 
thus valid for all the operators considered here even when n = oo. 

Unfortunately, H has almost never a purely discrete spectrum. It is well-known that 
the spectrum in general has continuous parts and resonances embedded in these, and a 
proper theory should treat these cases as well as the discrete part. In fact, the treatments 
of ifeff in the literature invariably glosses over this. It is an interesting future project to 
develop a geometric theory for the effective Hamiltonians which incorporates resonances 
and continuous spectra. 

The geometrical view simplified and unified the available treatments in the literature 
somewhat, and offered further insights into the effective Hamiltonians. Moreover, the the 
symmetry considerations in Sec. IIII Gl may have significant bearing on the analysis of per- 
turbation expansions and the properties of sub-cluster approximations to H^^. 

Indeed, it is easy to see, that if we have a complete set of commuting observables (CSCO) 



271] for Hq, and the same set of observables form a CSCO for Hi, all eigenvalues and 



eigenfunctions of H{z) = Hq + zHi are analytic in 2; G C, implying that the Rayleigh- 
Schroedinger perturbation series for H = Hq + Hi converges (i.e., aX z = 1) 2^. Intuitively, 
the fewer commuting observables we are able to identify, the more likely it is that there are 
singularities in H^si^z), so called intruder states. The Rayleigh-Schroedinger series diverges 
outside the singularity closest to 2; = 128}] , and in nuclear systems this singularity is indeed 
likely to be close to z = 0. On the other hand, resummation of the series can be convergent 
and yield an analytic continuation of H^s outside the region of convergence 29|]. To the 
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author's knowledge, there is no systematic treatment of this phenomenon in the hterature. 
On the contrary, to be able to do such a resummation is sort of a "holy grail" of many-body 
perturbation theory. A geometric study of the present kind to many-body perturbation 
theory and diagram expansions may yield a step closer to this goal, as we have clearly 
identified the impact of commuting observables on the principal vectors of £ and V. 

We have also discussed a compact algorithm in terms of matrices to compute ifgg, relying 
on the SVD. To the author's knowledge, this algorithm is previously unpublished. Since 
robust and fast SVD implementations are readily available, e.g., in the lapack library, 
and since few other matrix manipulations are needed, it should be preferred in computer 
implement at ions . 

As stressed in the Introduction, the algorithms presented are really only useful if we com- 
pute the exact effective Hamiltonian, as opposed to a many-body perturbation theoretical 
calculation, and if we know what exact eigenpairs to use, such as in a sub-cluster approx- 
imation. In this case, one should analyze the error in the approximation, i.e., the error in 
neglecting the many-body correlations in Hl^. In the perturbative regime, some results exist 
61] . The author believes, that the geometric description may facilitate a deeper analysis, and 
this is an interesting idea for future work. 
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